oo 

o 

O 
(N 

< 

en 
(N 



AC magnetization transport and power absorption in non-itinerant spin chains 
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We investigate the ac transport of magnetization in non-itinerant quantum systems such as spin 
chains described by the XXZ Hamiltonian. Using linear response theory, we calculate the ac mag- 
netization current and the power absorption of such magnetic systems. Remarkably, the difference 
in the exchange interaction of the spin chain itself and the bulk magnets (i.e. the magnetization 
reservoirs), to which the spin chain is coupled, strongly influences the absorbed power of the system. 
This feature can be used in future spintronic devices to control power dissipation. Our analysis 
allows to make quantitative predictions about the power absorption and we show that magnetic 
systems are superior to their electronic counter parts. 

PACS numbers: 75.10.Jm, 75.40.Gb 
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Power dissipation is one of the most important limita- 
tions of state-of-the-art electronic systems. The same 
is true for spintronic devices in which spin transport is 
accompanied by charge transport. In non- itinerant quan- 
tum systems, the dissipation problem is reduced since 
true magnetization transport generates typically much 
less power than charge currents [l|, |3] . This is one of the 
main reasons for putting so much hope and effort into 
spin-based devices for future applications [3|, y, |5|] . 

Here, we analyze non-itinerant quantum systems de- 
scribed by a spin Hamiltonian in which ac magnetization 
transport occurs via magnons or spinous (without the 
transport of charge). In Ref. [6[, the spin conductance 
of such a device has been derived with a particular fo- 
cus on the role of the magnetization reservoirs to which 
a one-dimensional spin chain is attached. We generalize 
this theory to the response to an ac magnetization source. 
This allows us to directly calculate (and thus estimate) 
the power absorption of such magnetic systems at a given 
driving frequency to using linear response theory. In gen- 
eral, the exchange coupling J in the spin chain and in 
the reservoirs will be different which is schematically il- 
lustrated in Fig. [TJ It turns out that the difference of 
the exchange coupling plays a crucial role in the depen- 
dence of the absorbed power as a function of uj. The 
larger the difference the stronger will be the suppression 
of power dissipation at finite frequencies. At low frequen- 
cies, however, the dissipative power is independent of the 
difference of the exchange couplings and takes a universal 
value determined by J in the reservoirs. 

We analyze the ac transport problem in quantum spin 
chains by a mapping of the spin Hamiltonian coupled to 
magnetization reservoirs to the so-called inhomogeneous 
Luttinger liquid (LL) Hamiltonian [7|,|8|,|9|. Interestingly, 
the absorbed power that is derived in this letter has an as- 
tonishingly simple dependence on the interaction param- 
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FIG. 1: (color online) Schematic of a quantum spin chain cou- 
pled to magnetization reservoirs. The magnetic field bias AS 
changes periodically in time. In the upper part of the figure, 
we illustrate that the exchange coupling in the spin chain Jsc 
(for la::j < L/2) can in general be different from the exchange 
coupling in the grey-shaded magnetization reservoirs Jrcs- As 
suggested in Ref. [y], this setup can be realized by a bulk ma- 
terial with an intra-chain exchange much stronger than the 
inter-chain exchange, where the material is heated to a tem- 
perature T > Tjv in the central part and cooled to T ^ Tn 
in the reservoir parts. (Tjv is the Neel ordering temperature.) 



eters of the LL model, see Eq. (|11|) below. This makes 
it a prime candidate for the experimental observation of 
LL physics in nature. 

In order to describe the system shown in Fig. [1] we con- 
sider a one-dimensional XXZ spin chain in the presence of 
a time-dependent magnetic field B(xi, t) — Bi{t)ez which 
can be described by the Hamiltonian H = i/xxz + Hb (t) , 
where 



H 



XXZ = J/l \^hxSj,x + Si,ySj^y + Asi^zSj.z) , (1) 



{^J> 



Hsit) = getJ'B^Bi{t)si^z. 



(2) 



Here, Si^a is the a-component of the spin operator at x^, 
(i, j) denotes nearest neighbor sites, ge is the g-factor, /is 



Bohr's magneton, and we assume anti-ferromagnetic cou- 
pling with J, A > 0. A possible realization of spin chains 
described by i?xxz is, for instance, a bulk structure of 
KCuFa or Sr2Cu03, where the exchange among different 
chains in the crystal is much weaker than the intra-chain 
exchange 
nian Hyrx 
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ly, lUl ll2 1 ■ It is well-known that the Hamilto- 
can be mapped onto a LL of spinless fermions 



Hll = -^ I dx 



9 



(3) 



where we have ignored Umklapp scattering [id] and made 
the identifications v — vs/g, vb — Jasm{kBa)/h, and 
g = (l + 4A/7r)"i/2 (a is the lattice constant). In Eq. ^, 
(p{x) is the standard Bose field operator in bosonization 
associated with spinon excitations here, n(x) its conju- 
gate momentum density, v the spinon velocity, vb the 
bare spinon velocity (at A = 0), ks the bare spinon 
wave vector, and g the interaction parameter (5 = 1 cor- 
responding to a non-interacting system, i.e. A = 0, and, 
in general for a -ffxxz spin chain, 1/2 < g < 1) ITlllSl]. 
In order to be able to properly describe the effect of reser- 
voirs, we modify the Hamiltonian i^LL in the spirit of the 
inhomogcncous LL model [7|,|8|,|9| described by a Hamilto- 
nian -ffiLL, where we assign a spatial dependence to v and 
g such that v{x) = vi and g{x) — gi being the spinon ve- 
locity and the interaction parameter in the reservoirs (for 
|a;| > L/2), respectively, and v{x) = v^ and 17(2;) = g^ 
being the corresponding quantities in the spin chain re- 
gion (for |a;| < L/2). Within this model, non-equilibrium 
transport phenomena such as the non-linear I — V char- 
acteristics and the current noise in the presence of irnpu- 
rities have been analyzed extensively 19|, |20|, |2l|, [23, l23| . 
In this letter, we are interested in a different situation, 
namely the ac magnetization transport in the linear re- 
sponse regime which should be seen coniplementary to 
the electric ac response analyzed in Ref. [2J, |25| . 
The Hamiltonian HB{t) describes a spatially varying 
and time-dependent magnetic field 5B{x) cos{ujt)ez with 
5B{x) = -AB/2 (AB/2) for x < -L/2 {x > L/2). For 
|a;| < L/2, 5B{x) interpolates smoothly between the val- 
ues ±AB/2 in the reservoirs [26|. The dc (w = 0) mag- 
netization transport of such a system has been analyzed 
in Ref. [6] and a spin conductance Gs = gi{ge^iBy' /h has 
been predicted. 

The magnetization current in linear response to an oscil- 
lating magnetic field can be evaluated using the following 
expression 



-'m v*^, ^j 



with 



dT 



dyaoix,y,T)dy6B{y,t-T) (4) 



ao{x, y, t) = 2i^^^^Q{r)dr ( Mx, T),^{y, 0)] ) (5) 



and the expectation value is taken with respect to -ffiLL- 
For X > L/2 and \y\ < L/2, the spin conductivity is given 

by 
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where ^ = {gi — gw)/ {gi + ffto) is the reflection coefficient 
of spinon excitations at a sharp boundary with different 
interaction coefficients gi and gw (Sl and 6(r) the Heav- 
iside function. The resulting spin current under continu- 
ous wave radiation reads 



4™)(x,t) = 2(l-7) 



givw [ge^-BY AS 



E^^^'x 



sin| )< cos 



UJ t 



p=0 

{2p+l/2)L L-2x 

Vw 2vi 

{2p + 3/2)L L-2x" 



LU t 



2vi 



0]} 



(7) 



We clearly observe an interaction-dependence of the mag- 
netization current in Eq. ([7]) through gi and 7. The pres- 
ence of higher harmonics due to higher order terms in 
7™ would be a strong experimental evidence for the spa- 
tial inhomogeneity of spin-spin coupling in realizations 
of XXZ spin chains. The physics behind the result in 
Eq. ([7]) is the following one: the system is driven with 
a continuous wave due to the ac magnetization source; 
therefore spinon excitations constantly enter and leave 
the spin chain from and to the reservoirs. Whenever, 
they experience a boundary in the exchange interaction, 
they are partly transmitted and partly reflected with a 
reflection coefficient 7. The resulting expression ([7]) is 
the superposition of all possible contributions to the spin 
current after infinitely many reflection processes. 

As a natural consequence, one may wonder whether 
an initial magnetization signal is actually transmitted 
through the spin chain. This depends crucially on the 
value of 7. To answer this question, we look at the mag- 
netization current in linear response to a unit pulse de- 
scribed by dy5B{y,t) = SBpSTpS{t ~ tQ)S{y — j/o) with 
2/0 G [—L/2, L/2] (where 6Bp corresponds to the height 
and 5Tp to the duration of the pulse). If we plug this 



expression into Eq. Q, we obtain for the spin current 



I^-'Hx,t) ^ gi6Bp6T, 



idefJ-Bf 



puip 



(i-7)e(i-to)E 



p=0 



a=± 



x-L/2 ^ L/2 - yo ^ 2pL i 



^W^sU + a 

x-L/2 L/2 + yo {2p+l)L 
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vi 



)}■ 

(8) 



The derivation of Im'{x,t) demonstrates that the ini- 
tially sharp 5-pulse is decomposed into a sum of in- 
finitely many J-pulses. Importantly, the amplitude of 
these pulses decreases by a factor 7 in a stepwise fash- 
ion once in each time interval L/v^ corresponding to the 
transit time in the wire. So, to answer the question how 
much signal has been transmitted we have to fix x, j/Qj 
to, and 7 in Eq. ([8]) and sum up all the prefactors of the 
(5-functions that can be non-zero in a given time inter- 
val between to and t. This analysis implies that all the 
dissipation happens in the leads and intrinsic relaxation 
is absent which is a direct consequence of the fact that 
the LL Hamiltonian describes a free boson. Once we in- 
troduce impurities the situation is different and intrinsic 
dissipation matters which we will briefly address below. 
We now turn to the discussion of the power absorption 
under continuous wave radiation. It turns out that the 
absorbed power is an ideal physical quantity to measure 
the reflection coefficient 7. We derive the absorbed power 
of the ID spin chain using Fermi's golden rule and linear 
response theory where particular care has to be taken 
of the spatial inhomogeneity of systems. The resulting 
expression is 



1 f /-^/^ 



L/2 



L/2 



dyY{.eao{x,y,uj) 



AB 



(9) 



where 

Re (Tq{x, y,uj) = ffu, " ^ |cos(a)(:r - y)) (10) 



h 

27(1 — 7^) cos((I') cos((I'(5: 4- y)) 
^ 1 + 74 - 272 cos(2w) 

27^ cos{i~j{x — y))(cos(2a)) - 7^) -1 
^ 1 + 74 - 272 cos(2(I') J' 

and we have introduced dimensionless variables x = x/ L, 
y = y/L, and uj = uj/uil with w^ = v^/ L. It is straight- 
forward to do the two remaining integrals in Eq. @ and 
the final result reads 



W{l 



9' 



{geliB^Bf /sin(c:;/2 



2h V '^/2 

l-74-H27(l-72)cos((I') 
1 + 74 - 272 cos{2Cj) 



(11) 




FIG. 2: (color online) The absorbed power is shown in units 
of Wo = (ge/is A_B)^/2/i as a function of the ac frequency u> in 
units of the ballistic frequency lul = Vw/ L. It is clearly visible 
that stronger repulsive interactions inside the wire decrease 
the absorbed power in the system. 



This is the main result of our work. It demonstrates that 
a measurement of the absorbed power due to ac response 
of the quantum spin chain is a simple and feasible way 
to measure interaction dependent coefficients such as g^ 
and 7. In Eq. (fTTj) . these coefficients just appear as pre- 
factors and not in a complicated power-law fashion as it 
is usually the case in observable quantities of systems de- 
scribed by LL physics. In Fig. [21 we show the interaction 
dependence of the absorbed power W{ijj). This demon- 
strate that stronger repulsive interactions inside the wire 
with respect to the leads suppress the dissipative power. 
If we compare Eqs. (O and pT|) we observe an inter- 
esting finite-size effect, namely that W{uj) vanishes as 
(sin(a)/2))2 close to a) « 27r whereas the leading contri- 
bution to Im {x,t) vanishes only as sin(a)/2) close to 
that driving frequency. Thus, the power absorption is 
more strongly suppressed than the magnetization current 
at frequencies close to 2ttll!l. This feature can be used 
in future devices to transfer data at special frequencies 
with low power dissipation. In the limit cj ^ 0, we obtain 
W{0) = gi{geHBA.B)'^ /2h corresponding to Joule heating 
where W{0) = I^{uj = 0)/(2G^). 

We now address the robustness of our main result pTjl 
against impurity scattering. An impurity can be mod- 
eled as an altered link in the -ffxxz chain, i.e. a local 
change in J on a nearest-neighbor link 271 128| . Within 
bosonization, such a scatterer at position xq in the sys- 
tem can be written as Hi = Xcos[\/ATT(p{xo,t) + 2kBXo]. 
If one of the two energy scales hvw/L or huj is larger than 
A (the bare impurity strength) , we can treat Hi perturba- 
tively up to lowest non-trivial order (which is second or- 
der in A) . In the presence of impurity scattering, the spin 
conductivity that enters into the calculation of Eq. © is 
subject to a (small) correction cr/(a;, y, lo) which has been 



derived for the corresponding electric case in Ref. [22| . 
For finite frequencies, ai{x,y,ijj) needs to be evaluated 
numerically. In the zero frequency limit, one finds power- 
law corrections to the spin conductance 29|, |30| resulting 
in power-law corrections to the absorbed power. In any 
case, as long as either hv^/L or huj are larger than the 
local change in J of the sample, the effect of impurity 
scattering is weak. 

The system which we considered previously consists of a 
spin chain smoothly connected to reservoirs. One may 
wonder how the previous result gets modified for isolated 
finite size spin chains, to which a time-dependent oscillat- 
ing magnetic field is applied along the chain (such that 
dB{x,t)/dx = A_Bcos(a;t)). For long Heisenberg chains, 
Hxxz still maps onto a LL of spinless fermions as in Eq. 
^ but with open boundary conditions (OBC). Follow- 
ing Ref. [3l|, we can estabhsh that Rea^^'^ {x,y,uj) = 
2ffm ^g'^^^ sm{ijjx / Vyj) sm{u!y/v^) for ui — ujn = TTnv^,/L 
{n — 1,- ■ ■ , {L~a)/a) and otherwise. From this expres- 
sion, we can infer directly the power needed to spatially 
shake the spin chain, using Eq. Q, 



W^OBc(^) =5t 



{geliB^Bf /sin((I;/2) 



sm' 



'('^/S), 



Cj/2 

(12) 
for (D = ujn/iOL and otherwise. Note that this power 
cannot be identified as dissipative power because a discon- 
nected LL does not contain a dissipative term. Instead, 
W^OBc('^) is the work per unit time needed to shake the 
system. This is the major difference to the case with 
leads, i.e. Eq. (fTTj) . where dissipation happens in the 
reservoirs. In the hmit a; ^ 0, IVqbc -^ due to the 
absence of reservoirs. 

Let us now compare typical values for the absorbed power 
in electric systems versus magnetic systems. We set 
9i ~ 9w = 1 for simplicity but keep in mind how finite 
interactions change the power absorption according to 
Eq. (jlip . The absorbed electric power in the dc limit is 
given by Wc\ — {e/WY /h. For a typical electric bias of 
Al/ = ImV, we obtain Wei ~ 3.87 • 10~"Js"i whereas 
the absorbed magnetic power for a typical magnetic bias 
of AS = O.IT is W « 2.59 • lO-^^Js"! (assuming g^, = 2) 
which is four orders of magnitude smaller. The rule of 
the thumb is Wci(Ay = O.lmV) ~ W(AB = IT). Thus, 
we expect substantial advantages of magnetic systems 
versus electric systems as far as power consumption is 
concerned. 

In summary, we have analyzed the magnetization current 
and the power absorption of quantum spin chains coupled 
to magnetization reservoirs with a time-dependent mag- 
netic field applied to the reservoirs. Both physical quan- 
tities depend crucially on the difference of the exchange 
interactions within the wire as compared to the magneti- 
zation leads. In fact, we envision to use this dependence 
as a way to control power dissipation in non-itinerant 
quantum systems in which magnetization transport oc- 



curs via spinous. Finally, we have briefly described the 
case of a finite size chain and the influence of impurity 
scattering on spin current. 
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